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Abstract 

We give some representation about recent achievements in analysis of the M/G/l queue 
with egalitarian processor sharing discipline (EPS). The new formulas are derived for the 
j-th moments (j £ N) of the (conditional) stationary sojourn time in the M/G/l —EPS 
queue with K (K £ 0 U N) permanent jobs of infinite size. We discuss also how to simplify 
the computations of the moments. 
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1 Introduction 

Processor sharing queues, made very attractive models by the works of Kleinrock 0, 0 and 
Yashkov 0.0. play a central role in queueing theory. These models were originally proposed to 
analyze the performance of scheduling algorithms in time-sharing computer systems, and continue 
to find new applications which pose interesting mathematical problems. Over the past few years, 
the processor sharing paradigm has emerged as a powerful concept for modeling of Web servers, 
in particular, for evaluating the flow-level performance of end-to-end flow control mechanisms 
like Transmission Control Protocol (TCP) in Internet. 

The mathematical analysis of processor sharing queues has resulted in many insightful results. 
Yet, a number of challenging problems remains to be explored. The main goal of this paper is 
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to gain understanding the problem of the moments of the stationary sojourn time in the M/G/l 
queue with egalitarian processor sharing (EPS), and to derive the formulas for the j-th moments 
(j E N) of the (conditional) sojourn time in the M/G/l—EPS queue with K (K 6 OUN) 
permanent jobs of infinite size. Our results complement and develop the corresponding sections 
of the paper by Yashkova and Yashkov ]Sj. 

An idea of the EPS discipline 1 was introduced by Klcinrock jT] who studied only M/M/1 case 
as a limit of the round-robin queue. In particular, he first showed that the mean sojourn time 
conditioned on the initial job size (service requirement) of the job is linear function of the size 
of the job. For an overview of the literature on processor-sharing queueing systems we refer to 
Kleinrock |2] (1976), Koboyashi and Konheim (1977), Jaiswal |7] (1982), Yashkov [3] (1987), 
[Hj (1990) and Yashkova and Yashkov |3}J (2003). 

The exact determination of the stationary sojourn time distribution in the M/G/l—EPS 
queue was an open problem for a long time. After puzzling researchers for 15 years, Yashkov 
(1981), [3j (1983) found an analytic solution of this problem in terms of double Laplace tranforms 
(LT) (all details contains also his book |XUj (1989)). Schassberger Q (1984) provided another 
(completely different) approach to the exact solution by considering the EPS discipline as a limit 
of the round-robin model (in discrete time). Later similar solutions were also made by means 
of the variants of the methods from [3] and from EH (or their combinations). See, for example, 
van den Berg P21 (1990) or Whitt [13] (1998) (here we do have no possibility to discuss the 
contributions of other authors (Brandt and Brandt TTj (1998), Asare and Foster (1983), Nunez- 
Queija (2000), Cheung et al. T3|(2005), et alii)to the closely related problems). We only mention 
that the EPS queue with permanent jobs has been studied in Hama Emu in from point of view 
which is different from our approach. A telecommunication system with CPU scheduling under 
SCO-UNIX can be considered as an example of using of the EPS model with permanent jobs for 
its description and predicting delays of the jobs. 

In fact, our method has turned out to be a very fruitful to derive many further results, for 
example, the time-dependent queue-length and sojourn time distributions in this and related 
models (see, for example, uni Ennuis] 2 ). These results hold for any stability condition. Besides, 

Wnder the EPS discipline, the processor (server) is shared equally by all jobs in the system. To put more 
concretely, when 1 < n < oo jobs are present in the system, each job receives service at rate 1/n. In other words, 
all these jobs receive 1/n times the rate of service which a solitary job in the processor would receive. Jumps of 
the service rate occur at the instants of arrivals and departures from the system. Therefore, the rate of service 
received by a specific job fluctuates with time and, importantly, its sojourn time depends not only on the jobs in 
the processor at its time of arrival there, but also on subsequent arrivals shorter of which can overtake a specific 
job. This makes the EPS system intrinsically harder to analyze than, say, the First Come — First Served (FCFS) 
queue. 

2 Indeed, the assumptions which required to use the steady-state (stationary) solutions of any queueing systems 
are rarely satisfied in real life. To be able really to apply queueing results in design and analysis of technical 
systems, in very many cases, the obtained results of steady state analysis are not sufficient. For example, it is 
often necessary to investigate the behaviour of the queue while it progresses towards a steady state (if and when a 
steady state exists). Even the average queue length at time t gives us much more information in comparison with 
the stationary mean of the number of jobs. 

However, few stochastic systems are known to have exact time-dependent (transient) solutions for the distri¬ 
butions of the processes. As a rule, such systems are the M/G/l queues with simpler disciplines (for example, 
FCFS, see, for example, Takacs jT7j). Besides, all time-dependent solutions of the queues of the type M/G/l 
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the entire transient and equilibrium behaviour of the M/G/l—EPS queue is contained in the 
results mentioned, and the most (if not all) available at present analytic solutions (and also many 
new) can be derived from them as special cases via standard arguments (for example, by means 
of the Abel’s/Tauber’s theorems). However, we shall not consider the transient solutions in this 
paper. 

The rest of the paper is organized as follows. In Section 2 we introduce some notations and 
describe our starting point represented by Theorem 12.11 In Section 3 we obtain some interesting 
consequences of Theorem 12.11 some of which were proved earlier as self-contained theorems but 
now ones are derived as special cases. The final section contains few closing remarks. 


2 Preliminaries 

In this section we give a short review of the M/G/l—EPS queue with K permanent jobs (only 
in steady state). For the time-dependent results we refer to [5] and also to EDI 13- 

jobs arrive to the single processor (server) according to a Poisson process with the rate A > 0. 
Their sizes (reguired service times) are i.i.d. random variables with a general distribution function 
B(x ) ((-£>(0) = 0, B( oo) = 1)) with the mean < oo and the Laplace-Stieltjes transform (LST) 
(3(s). Let f3j denote the j-th moment of B(x), j G N. The service discipline is the EPS: every job 
is being served with rate 1/n, when n > 0 jobs are present in the system. The EPS discipline is 
modified by having K > 0 extra permanent jobs with infinite sizes. The system works in steady 
state. In other words, p = \f3\ < 1 and very long time went from the instant 0 that marks the 
start of the work of our system till current time. 

ft is well known, due to Sakata et al. El that the stationary distribution (P n ) n >o °f the 
number of ordinary jobs in the M/G/l—EPS queue as K = 0 is geometrically distributed 

P 0 n = (l-p)p n , neOUN, (2.1) 

where p = A J 0 °°( 1 — B(x))dx < 1. We note that (Po n )n>o depends on the service time only 
through its mean. 

For K > 0 the equality m takes the form 

P Kn = (l-p) K+1 ( n ^ K )p n > neOUN. (2.2) 

We shall let that Vk(u ) denotes the conditional sojourn time of a job of the size u upon its 
arrival. This job enter into the EPS system with K > 0 permanent jobs in steady state. Let 
VKj = E[Vft'(u) J ]. (We shall omit the index K in these and similar notations when K = 0.) 

are obtained in terms of double transforms (on space and time) from which it is very hard to extract necessary 
information concerning the behaviour of the system. (Moreover, much more advanced mathematical techniques 
become necessary for the time-dependent solutions in comparison with steady state analysis.) Some exceptions 
give variations of the M/M/1—FCFS queue for which closed-form transient solutions are known. As a rule, the 
exact transient analysis of the M/M/1 —FCFS queue involves infinite sums of Bessel functions. In general, explicit 
exact solutions are highly unlikely for the time-dependent cases. 
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Define the LST of Vk(u ) by vk ( j , u ) = E[e~ rV ’ A ' ( “)j for Re r > 0 and u > 0. Let 7r(r) be the 
LST of the busy period distribution (due to ordinary, that is, non-permanent jobs). In other 
words, it is the positive root of the well-known Takacs functional equation IQ 


7r(r) = j3{r + A — A7r(r)) (2.3) 

with the smallest absolutely value. 

It is known from |2J the following theorem 


Theorem 2.1 When p < 1, 



vk(t,u) = E[e ryjf ( u )] = v(r,u) K+1 , 

(2.4) 

where v(r , 

u ) is given by the the equality \2.5l): 



v(r,n) = y 

ipfr, u) — a(r, 0, u) 

(2.5) 

Here 

poo 

a(r, 0,w) = A f>(r,u) * [e _ “^ r+A - ) (l — B(u))] + \e~ u( - r+ ^ / (1 — B(x))dx , 

(2.6) 

where is the Stieltjes convolution sign (on variable u), and fi>(r,u) is the LST (with respect 

to x) of some function ^(x,u) of two variables (possessing the probability density on variable x), 
which, in turn, has a Laplace transform (LT) with respect to u(argument q) 


it \ q + r + A f3(q + r + A) , 

f{r,q)= r>0,q> Avr r . 

(q + r + X)(q + A (3{q + r + A)) 

(2.7) 

In \2. 71), fi(r) = J“e rx dB(x) and n (r) (in the conditions imposed on \2.1\)) is understood as 
the minimal solution of the functional equation \2.3 1). 


Thus, the function f>{r, q) is given in the form of the two-dimensional transform of the function 
T(a;, u )) 


ip{r,q) = 


~ qu d x fy(x, u)du. 


( 2 . 8 ) 


'0 Jo 


In other words, ip(r, u) in equality (12.51) is the Laplace transform inversion operator, f>(r, u) = 
C^ l (f{r,q)){r,u), that is, the contour Bromvich integral 


^+2OO+0 


f>(r, u) = -—: / 'i/j(r,q)e qu dq. 


2m 


' —zoo+0 


Remark 2.1 Briefly, we have derived the expression for E[e -rV " K ^] by writing the sojourn time as 
some generalized functional on a branching process (like the processes by Crump-Mode-Jagers) 
by means of simple extensions of (non-trivial) arguments from (9 1 0 /. Using the structure of 
the branching process, we found and solved a system of partial differential equations (of the first 
order) determining the components of a (non-trivial, too) decomposition o/V^(m). It leads to 
E[e _r Df(D] ( see a i so Remark \fn$) . 
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3 Results 


We showed in the Section 2 that the determination of the steady-state sojourn time distribution 
in the queue M/G/l—EPS with K permanent jobs is simple extension of the results from |3], 
||. However, the solution contains the Bromwich countour integrals. First we consider the case 
K = 0. Equivalent form of (12.511 (without contour integrals) is given in the following theorem. 


Theorem 3.1 Equivalent form of \2. ,51) (without the Bromwich countour integrals) is given by 


1 

v(r, u ) 


_ rr> 

5^— fn(«), 

n\ 


n =0 


(3.1) 


where 

£o(u) = l, Uu) = j^-^u n - 1 *W^ n ~ 1) *(u), n = 1,2,... (3.2) 

Here pW n_1 )*('u) is (n — l)-fold convolution of the steady-state waiting time distribution W{u) 
in the familiar M/G/l—FCFS system with itself (W°*(u) = 1 (u), W l *(u) = W(u)), the LST of 
W (■ u ) is given by the well-known Pollaczek-Khintchine formula as 


w(q) 


1 ~P 

i -pf{qV 


(3.3) 


where f{q ) = (1— (3(q))/(q/3i) is the LST of the excess of B(-), that is, F(x) = P 1 1 (1 — B{y))dy 

(F°*(x) = l(x), the Heaviside function, F u (x) = F{x)). 


Proof. We rewrite (Ea in the form of Theorem 3.2 from m (see also (5.5) in [Sj), namely 

(1 - p)5(r,u) 


v(r, u ) = 


1 - p5(r, u ) I/“ + (1 - F(u)) 


(Re r > 0), 


(3.4) 


where 

5(r, u ) = e _u ^ r+A ' ) //>(r, u ) (3.5) 

and F(x) is introduced in Theorem 13.11 To reach our aim, it is used the LT of 1 /S(r,u) with 
respect to u (argument q), which is found from (12.71) as if(r, q — r — A), r > 0, q > r + A —A7r(r) 
(cf. also the third line on p.8 in [2j). Now we obtain after simple algebra the following power 
series expansion of the LT of the function l/v(r,u), r > 0, u > 0 



e- qu —^—du = 
v(r, u) 



(3.6) 


where w(q) is given by (13.31) . We note that 


rw(q) 


(i -p)q 


<lasg>r + A — A7r(r), p < 1. Now it is 


easily to invert analytically (on argument q) each term of the power series in r (13.61) . The result 
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is given by EH whence it follows (EH). the right-hand side of which is the power series in r with 
coefficients f n {u)/n\. □ 


The idea of such approach goes back to Heaviside. Similar results are obtained in [00], [20] • 
In fact, the form of Var[H(w)] [0], [3| (see the equality (13.101) below) stimulates a guess about the 
possibility of such expansion. 

Remark 3.1 The formula for W n *(x) in can be represented in the following form 

w™(x) = (i- P )”f; ( k+n _~ 

k =0 k ,2 / 

R is done, for example, by inversion of w(q) n , where w(q) is given by \3.3l) . 

Remark 3.2 We note that the by-product of our analysis is the distribution function W(x) whose 
LST is given by m- However, the analysis of EPS queue gives the other quantity (corresponding 
to a non-probability measure) W°(x) = W(x)/( 1 — p). The form of the LST ofW°(x ) is well- 
known: w°(q ) = P n f n (Q)- Unlike W(x), W°(x ) is well defined for all p > 0 and x > 0. R 
can be shown that W°(x) < oo for all p > 0, x > 0 and for any £>(•) (despite on the fact that, for 
p > 1, W°(x ) —> oo as x —> oo). 

Theorem 3.2 Let v n [u) = E[H(n) n ], n = 1, 2,... Then it holds the following recursive formula 

v n {u ) = r nN Jn n _ i (n)^,(n)(-l)* +1 (3.7) 

i =1 


Proof. Because v{r,u ) is analytical function in r (in particular, in r = 0), we can use the 
Tailor series expansion of v(r,u ) for small r > 0 

2 3 

iy> iyt 

v[r,u) = 1 - YyUi(u) + y V 2 (u) - g|-u 3 (w) + • • • (3.8) 

The product of (13.81) and (13.11) gives 

2 

- yyM«) - &(«)] + t^M^) - 2vi(u)fi(u) + f 2 (u)] 

3 

- -[v 3 (u) - 3 v 2 ( u )£ 1 ( u ) + 3vi(u)f 2 (u) - f 3 (tt)] + ... = 0 

and it leads to EH after differentiating n times with respect to r and setting r = 0. □ 

In particular, the expressions for the first two moments of V ( u ) are: 

Vi(u) = E[V (n)] = u/( 1 — p) (3.9) 
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(this is well-known result due to Sakata et al. jTHl (1969)), 

2 

(! — P)“ Jo 


Var[V(w)] = v 2 {u) — vf(u) = 


- I (u — x)(l— W(x)) dx, 


(3.10) 


where W(x) is introduced in Theorem 13.11 and it is expressed as 

OO 

w(n = (i-p)5>"F”*M 


(3.11) 


n=0 


(other variables were introduced above). 

The formula (|3.21) implies that £i(w) = E[V(w)] in (13.91) . The formula for the conditional 
variance «H5D was hrst obtained by Yashkov [2]. The standard way for the computation of the 
moments is the following 

*(«> Y^-i")^’ " eR (312) 

However, the LST v(r,u) in Theorem rm is very hard to differentiate in r more than once (prac¬ 
tically almost impossible matter) since this LST has a rather complex form due to a highly 
complicated form of such constituents of Ei> as a and 0. Therefore Var[H(M)] is hrst obtained 
by solving an alternative system of differential equations (see, for example, m Chapter 2]) which 
are derived by analogy with the equations of H Section 2] or with the equations in the proof 
of J3. Theorem 4], These equations are simpler forms of equations from jlj. ]3| because ones are 
composed not for the LST v(r, u ) but only for the second and the hrst moments. Thus the formula 
(EZj) for the case n — 2 in Theorem 13.21 was derived 20 years earlier than the same formula for 
arbitrary integer n (see ra, BSD 

ft can be useful for asymptotic expansion of v n {u ) for small and large u in the spirit of 
such expansion for Var[H(u)] (see final section for some details). Such results for VarjV('u)] were 
obtained at hrst in [2] (1981) (see also (3J). 

Now we show how to extend the formulas (13.91) and (13.101) to the case when the M/G/l—EPS 
queue is modified by having K > 0 extra permanent jobs with infinite sizes. 

Theorem 3.3 In the setting above, 

nV K (u)] = ( ^ + 1) ^ , (3.13) 

Var [V K (u)\ = ^ ^ f (w - x)(l - W(x))dx, (3.14) 

(1 — P) Jo 

where W(u ) is the steady-state waiting time distribution in the M/G/l—FCFS queue, represented 
by the equality 1X771/ . 

Proof. In our case, the equality EH takes the form 

VKn{u) = lim(—l n ) — ll ^ r n ' U \ ne N. (3.15) 

r \0 Or n 
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The formula (13.131) follows directly from (Q by means of applying (13T5D as n — 1. 

Taking into account (13 .1 f)|) . the formula EH follows also from El by means of applying 
EH as n = 2 after some simple algebra. □ 

Remark 3.3 An alternative way to obtain \3.1j\) is the following. We can compose and solve the 
system of the partial differential equations (of the first order) which satisfy the second and the first 
moments of Vk{u). The variant of such equations is known from JSJ TO] as K = 0. We point out 
the following fact. These equations rely on a decomposition of the sojourn time of the (tagged)job 
with the size u that arrives to the EPS queue when n standard jobs are present with remaining 
service demands xi,...x n (a key ingredient of analysis). Denoting this conditional sojourn time 
by VKn(u; x\,... x n ), it holds 

n 

V Kn {u ; xi,...x n ) = (.K + 1 )D{u) + ^ <f>(ay, u), (3.16) 

2=1 

where all components are independent random variables. 

The random variable D(u ) constitutes a “main” component of the sojourn time: it has the 
distribution of the sojourn time of a job with the size u that enters into a empty (from the 
standard jobs) system. By the way, its LST is given by \3.5\) . When the system is not empty, the 
i-th standard job (among the jobs which are sharing the capacity of the processor together with 
permanent jobs), having remaining size Xi, “adds” a delay <f>(oy, w) = <f>(ay Au,u) to the new job’s 
sojourn time. Note that D{u ) = <f>(ay, u ) for ay > u. Then the same chain of arguments as in 
can be used to derive 


4 Conclusion 


Using Theorems 12.11 and 13.21 we can easily obtain all other moments of Vk(u) in M/G/l—EPS 
queue with K > 0 permanent jobs. However, the exact expression even for the variance of the 
sojourn time (see EH involves an integration term, making an exact computations difficult 
from practical point of view. The same holds for the third, fourth, etc. moments. This difficulty 
remains also in the case K = 0. To overcome the difficulty, it is possibly to obtain some simple 
approximations for the second moments, see, for example, Villcla et al. [2] or van den Berg 
H2|. We note that there exist also an upper (and lower) bounds for Var[U(w)]. These bounds 
only depends on p and the size of job u. In addition, the bounds have the attractive property of 
intensitivity to B(x), and the difference between the upper and lower bounds is small, particularly, 
for small and moderate values of p. These second moments tight bounds can be easily generalized 
into higher moments of V(u) and also to the case K > 0. 

It are also known the asymptotic estimates of Var[U(w)] as u —> 0 and u —> oo [2]2i. For 
example, 


Var[V(u)\ 


u 2 p 


as u 


0 . 


(1-P) 2 

This is some asymptotics of the sojourn time variance of a very small jobs, and it leads to 
intensitive upper bounds with special structure requiring only knowledge of the traffic load and 










the job size. Now our results may be easily extended to the higher moments and also to the case 
K permanent jobs. Moreover, some preliminary analysis of asymptotics (see J9j) tells us about a 
high accuracy of such estimates in many typical cases. 
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